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METHOD AND APPARATUS FOR BACKSCATTER SPECTROSCOPY 

The present invention relates to a method of determining of a physical feature of a 
S medium, comprising: 

- producing radiation with a light source; 

- placiag a probe on a sample of said medium, the probe comprisiag a first optical fiber 
haviag a first diameter, and at least a second optical fiber having a second diameter; 

- sending light coming from the light source, through the jQrst optical fiber; 

< 

10 - collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber; 

- producing a first signal based on the first backscattered radiation^ and a second signal 
based on the second backscattered radiation; 

- determine a measured differential backscatter signal as a function of wavelength using 
1 S the first and second signals. 

Such a method is known from Amelink et al [1]. There, a special device is used to 
determine particle sizes in superficial layers. The device is suitable for measuring 
particle sizes in for example an aqueous suspension with polystyrene spheres, but is not 
fitted to accurately measure particle sizes in living tissue. So, detenrdning whether 
20 living tissue is normal or precancerous, by way of measuring particle sizes in living 
tissue is not very promising. 

In Doombos et al [2] the optical properties of human tissue are detemiined in 
vivo using a spectroscopic arrangement with ten optical fibers. One of the fibers is used 
to irradiate a sample, and nine other fibers collect the reflected light. By using a 
25 multitude of fibers to collect the reflected light, it is possible to calculate scattering and 
absorption coefficients of the sample. However, the method is not suitable for locally 
measuring the optical properties of the tissue. In particular, only mean values of the 
absorption coefficient of a relatively large part of the sample can be determined. 

It is an object of the present invention to locally measure a physical feature, such 
30 as a concentration, of a substance in a medium. 

The object is achieved by a method as described above, characterized by 

- calculating the physical feature by curve fitting said measured differential backscatter 
signal to a backscatter function, in which the backscatter function is a function of an 
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average path-length travelled by detected scattered photons, wherein the average path- 
ler^th is independent fironx an absorption coefficient of the medium, and from a 
scattering coefficient of the medium. Contrary to methods using diffiisely scattered 
photons such as Doombos et al [2], in the method according to the invention, the local 

5 absorption coefficient of the sample is measured in an absolute way, independent of the 
magaitude of the local scattering and absorption coefficients. This facilitates the 
measurement of absolute concentrations of absorbing molecules in a sample without 
requiring prior knowledge of the magnitude of the scattering and absorption 
coefficients of the medium. 

10 In an embodiment, the average path-length is proportional to the iBrst fiber 

diameter. This has as additional advantage that the average path-length and thereby the 
average penetration depth into the sample of the photons that contribute to the 
differential backscatter signal can be controlled by choosing the fiber diameter. As a 
result, the sampling volume can be controlled by adjusting the fiber diameter. 

15 Hence, the fiberoptic probe can be engmeered to match the relevant dimensions of the 
medium under investigation. 

In a particular embodiment, the physical feature is a concentration of at least one 

substance in the medium. 

The invention also relates to a device for determining a physical feature of a 

20 medium, comprising: 

- a light source for producing radiation; 

- a probe with at least a first and a second optical fiber, the first optical fiber having a 
first diameter and being arranged to deliver the radiation on a sample of said medium 
and to collect first backscatfcered radiation ftom said sample, the second optical fiber 

25 having a second diameter and being arranged to collect second backscattered radiation, 
wherein the second optical fiber is positioned alongside the first optical fiber; 

- a spectrometer for producing a first signal based on the first backscattered radiation, 
and for producic^ a second signal based on the second backscattered radiation; 

- a processor arranged to determine a measured differential backscatter signal as a 
3 0 function of waveler^th using the first and second signals, 

characterized in that the processor is arranged to calculate the physical feature by curve 
fitting the measured differential backscatter signal to a backscatter fimction, in which 
the backscatter fimction is a fimction of an average path-length travelled by detected 
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scattered photons, the average path-length being independent from an absorption 
coefGcient of the medium, and from a scattering coefficient of the medium. 

Furthermore, the invention relates to a computer program according to claim 
Sand to a data carrier according to claim 9. 

In another aspect of the invention, the invention relates to a method of 
determining a physical feature of a medium, comprising: 

- producing radiation ^th a light source; 

- placing a probe on a sample of the medium, the probe comprising a first optical fiber 
having a first diameter, and at least a second optical fiber having a second diameter; 

- sending light coming from the light soxirce, through the first optical fiber; 

- collecting first backscattered radiation through the first optical fiber and second 
backscattered radiation through the second optical fiber; 

- producing a first signal based on the first backscattered radiation; and a second signal 
based on the second backscattered radiation; 

- determine a measured differential backscatter signal as a fimction of wavelength using 
the first and second signals, 

characterized by 

- calculating the physical feature by curve fitting the measured dijQferential backscatter 
signal to a backscatter function, in which the backscatter fimction is a function of a 
mean firee path of photons. In this method, it is assumed that only singly scattered 
photons contribute to the differential backscatter signal and as a result the backscatter 
fimction can be easUy derived analytically. 

In an embodiment, the physical feature is a concentration of at least one 
substance in the medium. 

The invention also relates to a device according to claim 13. 

Furthermore, the invention relates to a computer program according to claim 14 
and to a data carder according to claim IS. 

Finally, the invention relates to a method according to claim 16. 

The present invention will be described below with reference to exemplary 
embodiments and the accompanying schematic drawings, in which: 
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Fig. 1 is a schematic diagram of a measuring device according to a preferred 
embodiment; 

Fig. 2a and 2b show cross sections of a sample and two fiber tips in the situation 
wherein the mean free path of the photons is much larger than the diameter of the 
5 fibers; 

Fig. 3 shows the results of Monte Carlo simulations for a homogeneous medium; 
Fig. 4 a dififerential backscatter signal normalized at zero absorption for several 
scattering coefficients; 

Fig. 5 shows a differential backscatter signal of a dilute suspension of 0.2 pm 
10 polystyrene spheres 

Fig. 6 shows the total differential backscatter signal as a fimction of the reflection 
coefficient ^(l) in the range of 1 0-1 00 mm"^ ; 

Fig. 7 shows the measured and calculated average path length r as a fimction of the 
average scattering coefficient; 
IS Fig. 8 is a graph of measurement results for three different absorption coefficients /4z 
showing the dc-fiber signal I, the c-fiber signal J and the differential backscatter signal 
Rbs as a function of wavelength; 

Fig. 9 shows a typical spectrum of an absorption curve A along with the specific 
absorption coefficient of Evans Blue dye; 
20 Fig. 1 0 shows a measured -4 * as a function of the absorption coefficient at A = 600 

9 

Fig. 1 1 shows typical spectra measured in a suspension of 1 .0 |Lun polystyrene spheres 
with and without Evans Blue dye; 

* 

Fig. 12 graphically shows a molar extinction coefficient as a function of wavelength; 
25 Fig. 13 shows in vivo measurements and a fit of the differential backscatter signal Rts 
in a human trachea realized using a fiber diameter of 400 pm, and 
Fig. 14 shows in vivo measurements and a fit of the differential backscatter signal Rts 
in a human trachea showing very low oxygenation indicative for lung tumor. 

30 A schematic diagram of a preferred embodiment according to the invention is 

shown in figure 1. The setup consists of a set of optical fibers for the delivery and 
collection of light to and from a sample 1 under investigation. Light from a light source 
2, for example a Tungsten Halogen lamp (Avantes HL-2000-FHSA), is led through a 
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first ann 3 of a bifurcated optical fiber. The bifurcated optical fiber is at a distal end 4 
coupled to a first distal end of a delivery- and-collection fiber 5 (in the following 
referred to as dc-fiber 5) which is small enough to be fit through a working channel of a 
clinical endoscope, not shown. A second distal end of the delivery- and-collection fiber 
5 5 contacts the sample 1 . Alongside the dc-fiber 5, a collection fiber 6 is arranged to 
collect light reflected by the sample 1. The collection fiber 6 (referred to as c-fiber 6) is 
connected to a slave channel of a dual-channel spectrometer 7, for example an Avantes 
SD2000. Preferably, the dc-fiber 5 is polished at a small angle to reduce specular 
reflections. * 

10 Light reflected back firom the sample 1 into the c-fiber 6 is led directly into the 

slave channel of the dual-chaimel spectrometer 7. A second arm 8 of the bifurcated 
fiber is connected to a master channel of the dual-chaimel spectrometer 7. Light 
reflected into the dc-fiber 5 is coupled back into the bifiircated fiber, and reaches the 
dual-channel spectrometer 7 via the second arm 8 of the bifiarcated fiber. An output of 

15 the spectrometer 7 is connected to an input of a processor 9 which is arranged to 
analyze signals firom the spectrometer 7. 

If only the dc-fiber 5 would be used to deliver and collect light to and firom the 
sample 1, a large firaction of collected light is due to single backscattering firom small 
sample depths, see [1]. A single to multiple scattering ratio depends on the scattering 

20 coefficient and phase function of the sample 1 and on a diameter of the dc-fiber 5. The 
contribution of multiply scattered Ught to the signal of the dc-fiber 5 can be 
approximately determined by combining the signal of the dc-fiber 5 with a signal 
coming firom an additional fiber, i.e. Ihe c-fiber 6 mentioned above. 

In [4] a differential backscatter signal Rts as a fimction of the wavelength X is 

25 determined using a formula like 



V \^ White (^) ' I biack (^)) J white (^)'J biack ( ^) 



(1) 



J 



30 where is the signal from the dc-fiber 5 in contact with the sample 1, In(X) is the 

signal firom the dc-fiber 5 submersed in a fluid with an appropriate refi:active index (for 
tissue: water would be appropriate), /w/„/e (X) is the signal firom the dc-fiber 5 with the 
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probe-tip at a specific distance from a diffuse reflecting reference material with a large, 
preferably wavelength-independent reflectance coefficient (white spectralon) and/ifadt 
(X) is the signal from the dc-fiber 5 with the probe-tip at that same specific distance 
from a diffuse reflecting reference material with a small, preferably wavelength- 
5 independent reflectance coefficient (black spectralon). Furthermore, J (X) is the signal 
fi»m the c-fiber 6 in contact with the sample 1 and J^huen,iack Q) is the signal from the c- 
fiber 6 with the probe-tip at the previously mentioned specific distance from the 
white^lack spectralon. Finally, c is a calibration constant that depends on the distance 
between the probe-tip and the reference materials. 

1 0 According to the invention, the processor 9 is arranged to calculate the physical 

feature using a predefined mathematical model, the differential backscatter signal (Rbs) 
and a curve fitting mechanism. In an embodiment, the diameter of the fibers 5, 6 are 
selected depending on a mean free path (m^) of photons sent into the sample 1. It is 
noted that if the mean free path can not be estimated before selecting a fiber diameter, 

1 5 initially two arbitrary fiber diameters may be selected. Afler curve fitting the measuring 
results using two different mathematical models, it will show which model applies. 

Figures 2a and 2b depict fiber tips of the dc-fiber 5 and the c-fiber 6 in the 
situation wherein the mean fise palli (niQ)) of photons coming out of the dc-fiber 5, is 
much larger than a diameter dfiter of the fibers 5, 6. In an embodiment, the diameters of 

20 both fibers 5, 6 are of equal size, however it should be understood that other selections 
are possible. In figure 2a, lines 21 and 22 show an example of a path traveled by a 
detectable singly scattered photon. In figure 2b, lines 23, 24, 25 and Unes 23, 24, 26 
show two possible paths of detectable multiply scattered photons. All multiple 
scattering events occur at such large distances from the fiber tip of the fibers 5, 6, that 

25 the probability of detection of multiply scattered photons is roughly equal for the dc- 
fiber 5 and the c-fiber 6. The differential backscatter signal Ris(X) will now purely be 
determined by singly scattered photons. 



In an embodiment, the diameter of the fibers 5, 6 are selected so that mQ)>4}Ae^. 
30 In the predefined mathematical model of this embodiment, the differential backscatter 
signal Rhs(X) is an ejqjonential ftmction of two times the mean free path. Below, an 
e^lanation for this model is given. 
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In the absence of absorbers, the differential backscatter signal Rbs(^) is 
proportional to the local, superficial scattering coefficient fJs(^)=^Qsca(^)'? ' As : 



where Capp is an apparatus constant that depends anaongst others on the distance 
between the probe tip and the reference materials (black and white spectralonX/7(A,Q) 
is a function called the phase function where Q is the scattering angle, Qsca(Ji) the 
scattering efficiency, p the concentration of substances present in the sample 1, and As 
the area of a scattering particle. For example, using a fused silica fiber with nxmierical 
aperture NA= 0.22, the differential backscatter signal Rbs(^) can be approximated by 



where (p is the azimuthal angle and 6 is the polar angle. 

Figure 3 shows the differential backscatter signal RbsW of measurements (see 
dots) of a dilute suspendon of 0.2 ^rni polystyrene spheres along with a calculation (see 
curve 32) according to Eq. (3). In figure 3 Rbs(X) is shown using arbitrary units (a.u.). 
Also a value of Qradar where Qradar = 4n' p{X, \^OyQsca(^), is indicated in the figure. 
Figure 3 shows excellent agreement between the measurement (i.e. the dots) and the 
calculation, which indicates that if m!^>dfiber , the single scattering is indeed the 
dominant contributor to the differential backscatter signal RtsW as defined in Eq. (1). 

A singly scattered photon first travels from the tip of the dc-fiber 5 to a particle, 
and then (the same distance) back firom the particle to the tip of the dc-fiber 5 (or tip of 
c-fiber 6), see also figure 2a. So an average path length r(X) traveled by the measured 
single scattered photons is equal to two times the mean firee path niQ)(A,), i.e. 

t(X) = 2 • mSp(X) (4) 




(2) 



RbsW^ C 




- Jd<p- Jde.sin(e)./7(^i80).//,(A) ^c^p'piknoyMsa) (3) 



0 170 
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_ • 

In the presence of w absorbing species with specific absorption coefificients /Ja^^'^'(^)y 
the differential backscatter signal becomes 



RM = Capp' -paiSO) • /^A) • exp(- m • |] Pi - ^la'^'\X)) 



Ml 



n 



= Capp' • 180) • fi^X) ■ exp(-2 • m^(A)) • J] pi • nT'^XX)) (5) 



Ml 



where Capp' is an apparatus constant, p{X,l%G) is the phase function, /4(>l) is the 

scattering coefficient of the medium, X is the wavelength of the first and second 

backscattered radiation, mSgiX) is the mean free path as a function of the wavelength, n 
10 is the number of substances in the sample 1, pi is the concentration of absorber i present 

in a detection volume of the sample 1, and jLia^^^*\X) is the absorption coeflBcient of 

absorber / as a function of the wavelength. 

It is noted that in Eq. (S) the assumption is made that absorbers are 

homogeneously distributed and do not influence each other. The Eq. (5) may be 
1 S corrected for non-linear phenomena such as an inhomogeneous distribution of 

absorbers, see e.g. [8]. 

According to an embodiment, the specific absorption coefiicients of the 

absorbers, the wavelength dependency of the scattering coefficient jUs and the phase 

function /?, together with Eq. (5) are used in order to calculate the concentrations of all 

20 the absorbing substances present in the detection volume of the sample 1 . Since the 
detection volume is typically very small in the present invention, the extracted 
concentrations are highly spatially resolved. This is not possible with the known 
methods that are based on diffuse reflectance, and wherein the obtained concentrations 
are averages over large sample volumes, see e.g. [2], 

25 The apparatus constant Capp' (Eq. 3) can be detemiined for a specific distance 

between the tip of the dc-fiber S and the reference materials (black and white 
spectralon). For a suspension of monodisperse polystyrene spheres of known size and 
concentration, the scattering coefficient //^ and the phase function jp(180) can be 
calculated using Mie theory [4]. The apparatus constant Cg^^p' simply follows from Eq. 

30 (3). In terms of the volume fraction / of the suspension, the radius of the spheres a and 
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the radar efficiency coefficient QradmQ^'^ p(k ,18O)-0^cfl(^) the apparatus constant is 
determined by 

RbsQ^) « Capp • pQ^,\^0) MsQ^) = Capp • 0.05968 • - • QradafQ^) (6) 

a 

According to another embodiment, the selected diameter dflber is chosen so that 
the mean free path is smaller than dftber* In this embodiment, the differential backscatter 
signal Rbs is a function of the fiber diameter dfiber* This will be discussed in more detail 
below. 

When the mean free path of the photons is smaller than the selected fiber 
diameter (i,e. m£p(X)<dfibei)9 the contribution of multiply scattered light to the 
differential backscatter signal Rbs(^) of the single dc-fiber S cannot completely be 
removed using 6q. (1). In this case, it appears that the average path length of the 
photons contributing to the signal Rbs(^) becomes nearly independent of the optical 
properties of the sample 1. In this situation, multiple scattering events already occur at 
small distances from the tip of the dc-fiber 5. An analytical expression for the 
backscatter signal Rbs(X) is not available for this situation and Monte Carlo simulations 
were used to model the behaviour of RbsW as a function of the diameter of the fibers 5, 
6 and of the optical properties of the sample 1. Figure 4 shows the results of Monte 
Carlo simulations using the MCML-code (Monte Carlo for Multi-Layered media) of 
Wang et al [6,7] for a homogeneous medium with an anisotropy value g=0.9. A flat 
circular incident beam with diameter dflber is directed onto the sample 1, and the 
differential backscatter signal Rbs is calculated by subtracting the total reflectance in the 
c-fiber 6 (with diameter dflber and center located at a distance dflber fi^om the center of the 
incident beam) firom the total reflectance in the dc-fiber 5 (witli diameter dflber 
overlapping the incident beam). Simulations were performed for sets of four different 
scattering coefficients (m=15, 25, 50 and 80 mm"^), four different fiber diameters 
(^6er=200, 400, 600 and 800 ^mi) and five different absorption coefficients (|Lia=0, 0.2, 
0.4, 0.6 and 0.8 mm"^). 

Figure 4 shows Rbs as a function of absorption coefficient jia where the open 
circles/dashed lines correspond to ^6^^=200 pm, the filled circles/dotted lines 
correspond to ^^fer==400 fjm, the open squares/solid lines correspond to 4jfer=600 |xm 
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and the filled squares/dash-dotted lines correspond to dftber=^00 ]jin. The differential 
backscatter signal Rbs for each scattering coefficient \Xs was normalized to the 
(dftber^200 |j,m, jaa=0 mm'^) case. Figure 4 shows that in the absence of absorption, i.e. 
Ha = 0, the differential backscatter signal Ri,s depends linearly on the scattering 
5 coefficient jig. Furthermore, the slope of the straight lines (signifying the relation 
between Rbs and |ia) depends only on the fiber diameter and is independent of the 
scattering coefficient The latter is more clearly demonstrated in figure 5, where the 
differential backscatter signal Rbs is normalized to unity at zero absorption for all 
scattering coefficients i^s. The open circles correspond to dftber^2Q0 and the filled 
10 squares correspond to ^6er=800 |jm. These Monte Carlo simulations therefore suggest 
that in the situation where m^<dfiber> the diameter of the fibers 5, 6 determines the 
average path length x of the measured photons. The backscatter signal Rbs for this range 
of parameters can thus be written as 

15 Rbs(^) = Cr m • exp(-T • |J«) = Ci- i^s • e3cp(-C2 • dftber ' (7) 

where Ci and C2 are constants, x is the average path length, (Xa is the absorption 
coefficient, m is the scattering coefficient and dftber is the fiber diameter of the fibers 5, 
6. 

20 An exact analytical expression for Rbs(^) is not available due to the large 

contribution of multiple scattering events to the signal. Measurements were done for 
determining a total integrated backscatter signal Rtot for a range of A between 400-900 
nm, using tiie formula 

900Dm 

25 Rtot(Ps)= J dXRbsdM^) (8) 

400nni 

Figure 6 shows that the integrated total backscatter signal RtotM proportional to 
jUs(X) in the relevant range of between 10-100 mm"^ Therefore, in the absence of 
absorbers, it follows that 
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Rbs(^) = CappMs(X) (9) 

which is in agreement with the Monte Carlo simulations. 

5 

In the presence of n absorbing substances in a suspension, with specific 
absorption coefficients /fc^^'^Y'V » the differential backscatter signal becomes 

Rtsa) = C,5,y/<,(2;-exp(- r • Pi • Ma'^'^'a)) (1 0) 

10 

where r is the average path length of the detected backscattered photons and pi is the 
concentration of the substance i. 

Non-linear phenomena such as an inhomogeneous distribution of absorbers are 
not incorporated in Eq. (10), but can be added by the skilled person, see e.g. [8]. 

1 5 Figure 7 shows the measured and calculated average path length r as a function 

of the averse scattering coefficient <fis(X)> (with 500 nm < A, < 700 nm) for dfiher = 0.4 
mm and for absorption coefficient fJa(X) = 2.0 mm"^ at A=600 nm. In figure 7 
measurement are indicated by dots, the r = 2-mJ^ curve is indicated by a line 71 and 
Monte Carlo simulations are depicted by dashed lines. Identical results were obtained 

20 for suspensions with an absorption coefficient of = 1 .0 rmn'^ at 600 nm. The average 
path length r was determined using suspensions of polystyrene spheres with different 
sizes and concentrations to vary the scattering coefficient . The anisotropy g of 
these suspensions was in the range of 0.8-0.9. Evans Blue dye was added as an 
absorber, and the average path length t was calculated from Eqs. (9) and (10) and 

25 knowledge of the concentrations and specific absorption coefficient of Evans Blue, as 
will be known to the skilled person. 

Looking at the measxjred average path lengths of figure 7, it clearly shows that for 
large scattering coefficients (/ijr = 10-100 mm"\ the range relevant for tissue) the 
average path length r is independent of the scattering coefficient jl^ to within 10% and 

30 approxnnately equal to half the fiber diameter (r« 0.24 mm while dfiber = 0.40 mm), in 
agreement with the Monte Carlo simulations (the dashed lines correspond to Monte 
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Carlo calculations for dftber 0.2, 0.4, 0.6 and 0.8 mm). For small scattering coefficients 
(e.g. /4 < 5 mm"^) the average path length r is well described by t= 2-ni^, according 
to Eq.(4), see line 71. Figure 7 also clearly demonstrates that the transition from the 
'single scattering regime' to the ^constant path length regime' occurs for mean free 
5 paths of the order of the fiber diameter. It is therefore expected that single scattering 
will prevail over a larger range of scattering coefficients for fiber diameters smaller 
than 400 pm. 

In the following, the effect of absorption on the average path length r will been 
examined in more detail. Various concentrations of Evans Blue dye were added to a 

10 suspension of polystyrene spheres with scattering coefficients /4 of 35 mm"\ The 

concentrations of Evans Blue (EB) dye were varied such that the absorption coefficient 
fJa at 600 Em was in the range of 0 to 2 mm"^ Typical results of the differential 
backscatter signal R^s for three different absorption coefficients /Ja are shown in figure 
8. Note that the signal / of the dc-fiber 5 is plotted on a different vertical scale than the 

15 signal J of the c-fiber 6 and the differential backscatter signal i?^^. 

The spectra with Evans Blue present in the suspension were divided by the 
spectrum with no Evans Blue present R° and the negative natural logarithm of the ratio 
R™/R^ was determined: 



20 ^ = -hi(R™/R^ = x-p-/4,^^'™ (11) 



where p is the concentration of Evans Blue, and is the specific absorption 

coefficient of Evans Blue. 

25 Figure 9 shows a typical spectrum of an absorption curve 92, along with the 

specific absorption coefficient fj^^^"^^^ of Evans Blue dye, see curve 94. 

For all concentrations, an area A * under the absorption curve 92 was determined 
in the wavelength range % between 500 and 650 nm. From Eq. (1 1) it follows that if the 
average path length x is independent of the absorption coefficient /4,^^% area^* should 
30 depend linearly on the concentration p of the species in the suspension. 

Figure 1 0 ^ows a measured ^ * as a ftmction of the absorption coefficient Ha at 
600 nm. A curve fitting is done, resulting in a line 104 for ///A )=35 mm"^ Figure 10 
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shows that the average path length x is indeed independent of the absorption coefficient 
fia in the range 0-2 nun^^ 

From the previoxis results of figure 4 to 10, it shows that for < dfiher^ the 
differential backscatter signal i?^ is described by Eq. (7) with C2 « 0.6. Figure 1 1 
S shows typical spectra measured in a suspension of 1.0 |im polystyrene spheres with and 
without Evans Blue dye (|Xa= 2 and 0 mm"^ at 600 nm, respectively). From Eq. (7) and 
figure 7, Ihe relation between the differential backscatter signals with and without 
absorber is given by 



10 RhsiKHt^ =^ Rbs(^0yexp(-0.24 ^ia> (12) 



The calculated spectrum according to Eq. (12) is plotted as a dashed line 1 10 in 
figure 1 1 and shows excellent agreement with the measured RbJ^KMa)^ see line 1 1 1 . In 
figure 1 1 , line 1 12 depicts RbiJi^ 0). 
IS In short, the average path length t of photons measured when subtracting the 

signals of the c-fiber 6 from the dc-fiber 5 using Eq. (1), is independent of the optical 
properties of the sample 1 and approximately equal to half the diameter of the fibers 5, 
6 used, as long as the fiber diameter is larger than the mfp. 



20 In a specific embodiment, the device according to the invention is arranged to 

determine concentrations of oxygenated blood in tissue. Since the scattering coefficient 
of tissue ///^^^ is in the range of 10-100 mm^^ the fiber diameter should be smaller 
than a certain maximum diameter dmax where dmax is between 10 and 100 |am, so for 
example smaller than 50 mn, in order to measure predominantly single scattering in 

25 tissue. In this case, Eq. (5) holds. For fibers 5, 6 with much larger diameters (e.g. 200 
or 400 ^im), the differential backscatter signal RbsQ^) is described by Eq. (10) with r« 

It is presently known that the wavelength dependence of the scattering coefficient 
in tissue jLiJ'^ can be adequately described by an empirical power-law function, see 
30 also [3], [4], [5]. 



(13) 
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with a and b constaats that depend on the size, concentration and relative refractive 
index of the scatterers (i.e. substances) present in the detection volume. 
The dominant absorbers in tissue in the visible wavelength range are oxygenated and 
5 deoxygenated blood. Thus in tissue Eq. (10) becomes 

Rbiy^) = Capp • aA,-^ • exp(- O.&dflter • Pblood • (So2 • >t/«^^^^ + (I-S02) • )) 

= C'app • • exp(- O.&dj^er ' Pblood • (So2 * Ai^'^"^ + (1 -S02) ' juj^'^'^'' )) (14) 

10 

where pbiood is the concentration of blood, Sq2 is the blood oxygenation (percentage 
oxygen saturation) in a certain detection volume, Capp is a constant that depends on the 
caUbration constant c, Capp is Capp-a, X is the wavelength, b is the slope of the 
scattering coefficient defrtied in Eq. 13, fi^^^ is the specific absorption coefficients of 

15 fiilly oxygenated blood, ^jJp^^^^ is the specific absorption coefficients of fully 
deoxygenated blood. 

Non-linear phenomena such as an inhomogeneous distribution of absorbers are 
not incorporated in Eq. (14), but can be added by the skilled person, see e.g. [8]. 

Since the specific absorption coefficients of fiilly oxygenated (juj^^^^ and fully 

20 deoxygenated (/la^^^'^^^ blood are well known, see figure 12, Eq. (14) can be fitted to 
the measured data to yield the slope b of the scattering coefficient /z/^*^^, the 
concentration pbiood and the oxygen saturation S02 of the blood present in the detection 
volume. When a correction is made for the inhomogeneous distribution of blood in the 
vessels, a vessel diameter D may be determined as well. Since the average detection 

25 depth is small (e.g. 0. 1 nmi), fhe blood present in the detection volume when measuring 
non-invasively is located in capillaries. 

In figure 13, in vivo measurements of backscattering in a himian trachea together 
with a fit using Eq. (14) are shown. The measurements are realized using a fiber 
diameter of 400 pm. The dots depict the measurements and curve 130 is a fitting curve. 

30 In figure 13,6 = -0.94 and the oxygenation S02 = 95%. 

The present invention can be used for tumor detection. Tumor growth may, due 
to its excessive oxygen consumption, be accompanied by a low capillary oxygen 
saturation, which can only be assessed using a very localized measvurement. Since (pre- 
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)cancerous tissue is generally more heterogeneous than normal tissue, the standard 
deviation of multiple measurements is likely to be larger for (pre-)cancerous tissue than 
for normal tissue. Standard deviations in the measurements can be calculated for the 
oxygen saturation, the blood concentration, the blood vessel diameter, and the slope b 
5 of the scattering coefficient /z/"^^. It is noted that the invention is by no means 
restricted to determine a concentration of a substance as the physical feature. All 
features, mentioned in the previous phrase can be regarded as physical features. 

An example of a measurement of a Ixmg tumor is shown in figure 14. The shape 
10 of the dip in the wavelength range of 500-600 nm in this figure demonstrates the 
depletion of oxygen from the capillaries of this tumor due to its excessive oxygen 
consumption. 

When a needle-probe is used, the local oxygenation and scattering coefficient 
^f"^^ can be measured invasively. This could be helpfid in determining tumor-margins 

15 intra-operatively in real-tune, for instance during resection of a breast-tumor. 

According to an embodmient, the device comprises multiple probes and a 
multichannel spectrometer for multiple shnultaneous measurements on diSerent 
locations of the sample 1. Using this device, multiple measurements can be made 
simultaneously on differ^t locations of for example a suspicious lesion. 

20 In yet another embodiment, the device comprises at least two pairs of fibers, 

having different fiber diameters. For example, when a pair of fibers with 100 |xm, a pair 
of fibers with a diameter of 200 |jm and a pair of fibers with a diameter of 400 pm are 
used, information from different depths in the sample 1 can be obtained as the average 
path length increases with increasing fiber diameter. 

25 The method and apparatus according to the invention can also be used to analyze 

local dmg concentrations. From Eq. (10) it follows that if the specific absorption 
coefficient of a certain drug is known, the local concentration p of that substance can 
be determined usmg the invention. 

Another possibiUty of the present invention is to monitor glucose concentrations. 

30 The scattering coefficient /^^^^^ depends among others on tiie relative refractive index 
of the scatterers with respect to the surrounding medium (in tissue: cytoplasm). The 
refractive index of the surrounding cytoplasm is likely to depend on the concentration 
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of glucose. A chaage in the glucose concentration will therefore likely affect the slope 
b of the scattering coefficient pj^"^^ , see Eq. (13). 

Whilst specific embodiments of the invention have been described above, it will 
be appreciated that the invention may be practiced otherwise thaa as described. For 
5 example, a concentration of a substance in polluted water may be calculated. The 
description is not intended to limit the scope of the invention. 
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